Are Two z ~ 6 Quasars Gravitationally Lensed ? 
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Abstract 

^ ■ Several high-z (z > 5.7) quasars have been found during the course of Sloan Digital 

Sky Survey. The presence of such very high- z quasars is expected to put constraints 
^ ' on early structure formation. On one hand, it is suggested that these most luminous 

objects at high redshift are biased toward highly magnified objects by gravitational 
! lensing. To clarify the effect of gravitational lensing on high-z quasars, we began an 

S : imaging survey of intervening len s in g galaxies. Indeed, onr previous optical image 

showed that SDSSp J104433.04+012502.2 at 2=5.74 is gravitationally magnified by 
^ ■ a factor of 2. In this paper, we report our new optical imaging of two other high-z 

quasars, SDSSp J103027.10+052455.0 at z=6.28 and SDSSp J130608.26+035626.3 
at z=5.99. Since we found neither an intervening galaxy nor a counter image with 
%' < 25.4 — 25.8 around each quasar, we conclude that they are not strongly magnified 
regardless that a lens galaxy is dusty. 
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1. Introduction 



The Sloan Digital Sky Survey (SDSS: e.g., York et al. 2000) has been producing the 
discovery of high-redshift quasars at z ~ 6 (Fan et al. 2000, 2001, 2003); the most distant one 
known to date is SDSSp J114816.64+525150.2 at z = 6.43 (Fan et al. 2003). One significant 
feature seems to be that all of these z ~ 6 SDSS quasars are exceptionally bright. As claimed 
by Wyithe & Loeb (2002b; hereafter WL02), the presence of such very bright quasars at z ~ 6 
raises the following two problems: (1) If the central engine of these quasars is an accreting 
super-massive black hole, and shines near the Eddington accretion rate, a super-massive black 
hole with mass exceeding ~ 3 x 1O 9 M was already formed beyond z — 6. This challenges any 
theoretical models for early structure formation (Turner 1991; Haiman, Loeb 2001). (2) Also, 
the absolute luminosities of the z ~ 6 quasars are systematically higher than expected from the 
SDSS survey criteria (Fan et al. 2000, 2001), suggesting that the luminosity function of the 
z ~ 6 quasars is significantly biased to higher luminosities. These two problems led WL02 to 
propose that a third of the known z ~ 6 quasars may be magnified by a factor of 10 or more 
by intervening gravitational lenses, although the lensing probability depends on the slope of 
the luminosity function (Wyithe, Loeb 2002a). Comerford et al. (2002) also estimated the 
probability of gravitational lensing magnification by a factor of 10 or more for the four SDSS 
quasars at z ~ 6, to be ~ 100% given that the LF of the quasar is steep or extends to faint 
magnitudes. 

During the course of our optical deep-survey program on Lya emitters at z ~ 5.7 in 
the field surrounding the quasar SDSSp J104433.04-012502.2 at z = 5.74 (Ajiki et al. 2002; 
Taniguchi et al. 2003), we found that a faint galaxy with m£(AB) pa 25 is located at l."9 
southwest of this quasar (Shioya et al. 2002). We estimated the redshift of this galaxy to be 
z ~ 2.1, and concluded that SDSSp J104433. 04— 012502.2 could be magnified by a factor of two. 
Although this fact supports the suggestion by WL02, further information on other very high- z 
quasar is necessary. In this paper, we present our new results of a lensing galaxy survey around 
two SDSS quasars, SDSSp J103027.10+052455.0 at z = 6.28 and SDSSp J130608.26+035626.3 
at z = 5.99. Throughout this paper, we adopt a flat universe with Q m = 0.3, Q\ = 0.7, and 
h = 0.7, where h = H /(100 km s" 1 Mpc" 1 ) and we use magnitudes given in the AB system. 

2. Observations and Results 

2.1. Observations 

Deep optical imaging observations were made with the Suprime-Cam (Miyazaki et al. 
2002) on the 8.2 m Subaru Telescope (Kaifu 1998) at Mauna Kea Observatories. The Suprime- 
Cam consists of ten 2kx4k CCD chips and provides a very wide field of view: 34' x 27' with a 
0''2/pixel resolution. Using this facility, we made a very deep i'-band imaging survey in the fields 
around SDSSp J103027.10+052455.0 at a redshift of 6.28 and SDSSp J130608.26+035626.3 at 
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a redshift of 5.99 (Fan et al. 2001; Becker et al. 2001). A summary of the imaging observations 
is given in table 1. All of the observations were made under photometric conditions, and the 
seeing size was between l."2 and l."8 during the run. 

The CCD data were reduced and combined using IRAF and mosaic-CCD data-reduction 
software developed by Yagi et al. (2002). Since the observation was made during the course of 
the Subaru Deep Field project (e.g., Kodaira et al. 2003), the available observing time for the 
purpose was limited, and thus we could not take photometric-standard star data. Therefore, in 
order to calibrate our data, we used the SDSS Date Release 1 (DR1) catalog. We selected all 
point sources within 5 arcmin from each quasar from the SDSS DR1 catalog. In this procedure, 
we selected the point sources with 18 < %' < 20, because objects with %' < 18 are too bright in our 
images, while objects with i' > 20 are too faint to be well calibrated in the SDSS DR1 catalog. 

After the calibration, we obtained the %' magnitudes of each quasar: %' = 23.23±0.07 
for SDSSp J103027.10+052455.0, and i! = 21.95±0.03 for SDSSp J130608.26+035626.3. The 
magnitude of SDSSp J103027.10+052455.0 is consistent with that of Fan et al. (2001). On the 
other hand, the magnitude of SDSSp J130608. 26+035626. 3 is 0.63 mag brighter than that of 
Fan et al. (2001). The difference may be because our image of SDSSp J130608. 26+035626. 3 
has a deeper limiting magnitude, and thus our photometric accuracy is better than that of the 
SDSS data for faint objects. 

2.2. SDSSp J103027.10+ 052455.0 

The i'-band images around SDSSp J103027.10+052455.0 are shown in the left panel 
of figure 1. We did not detect any galaxy within 6 arcsec from the quasar SDSSp 
J103027. 10+052455.0. No possible counter image of the quasar was identified either. We 
checked the radial profile of the quasar in this image, and found that the image is consistent 
with that of point sources. Therefore, we conclude that the quasar may not be gravitationally 
magnified by a foreground galaxy. 

Fan et al. (2003) analyzed high-resolution images of SDSSp J103027. 10+052455.0 ob- 
served by the Advanced Camera for Surveys on the Hubble Space Telescope and the Near 
Infrared Camera on the Keck I telescope, and found that the quasar appears to be an unre- 
solved point source. They also mentioned that there are no close companions with K' < 21 
within 10" from the quasar. These results are consistent with our results. Schwartz (2002) 
pointed out that SDSSp J103027. 10+052455.0 is not consistent with a point source from the 
X-ray observation by the Chandra X-ray Observatory, and interpreted that the quasar may be 
gravitationally lensed. However, we found neither a galaxy nor a counter image within 6 arcsec 
from the quasar. This fact shows that the quasar is not strongly magnified. 
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2.3. SDSSp J 130608. 26+035626.3 

The i'-band image around SDSSp J130608. 26+035626. 3 is shown in the right panel of 
figure 1. Again, we found neither a galaxy within 3 arcsec from the quasar nor a counter image 
of the quasar within 5 arcsec around SDSSp J130608. 26+035626.3. We also checked the radial 
profile of the quasar in this image, and found that the image is consistent with that of point 
sources. 

Fan et al. (2003) also analyzed the high-resolution images of SDSSp 
J130608. 26+035626. 3 derived from the Astrophysical Research Consortium 3.5 m telescope, 
and concluded that the quasar also appears to be an unresolved point source. Ivanov (2002) 
obtained a deep /-band image of the quasar from the VLT with FORS1, and showed the radial 
profile of SDSSp J130608. 26+035626. 3 is consistent with that of a point source. Our result 
concluded that the quasar is an unresolved point source as well. Schwartz (2002) also pointed 
out that the quasar is consistent with a point source. Becker et al. (2001) found a candidate of 
damped Lya systems at z = 4.86 in the spectrum of the quasar. However, we found neither an 
optical counterpart nor a counter image of the system within 3 arcsec from the quasar in our 
i'-band image. We conclude that the gravitational magnification, if any, of the quasar is not 
strong. 

3. Discussion 

We cannot find any lens galaxy in our i' images around SDSSp J103027. 10+052455.0 at 
z=6.28 and SDSSp J130608. 26+035626. 3 at z=5.99. However, these results do not mean that 
the quasars are not magnified by gravitational lensing of intervening galaxies, since galaxies 
fainter than our limiting magnitudes may still be present and magnify the light of quasars. We, 
therefore, evaluated the maximum magnification factor by galaxies fainter than the limiting 
magnitude. In order to estimate how the quasar is gravitationally magnified, we need informa- 
tion about both the redshift and the stellar velocity dispersion of the lens galaxy. Although 
there is no information on z for a possible lens galaxy, we can estimate the maximum luminosity 
of galaxies which is consistent with our observations. Second, we also estimated the maximum 
velocity dispersion of the dark matter halo of galaxies using the Tully-Fisher relation. Then, 
we estimated the maximum magnification factor based on the singular isothermal sphere model 
for gravitational lensing. 

3.1. Absolute Magnitude of the Lens Galaxy 

The maximum luminosity, which is consistent with our observation, depends on the shape 
of the spectral energy distribution (SED) of a lens galaxy. The SED is mainly determined by 
the following four factors: (1) the radiation from stars in the galaxy, (2) the extinction by 
dust in the galaxy, itself, (3) the redshift of the galaxy, and (4) the absorption by intergalactic 
neutral hydrogen between the galaxy and us. We treat the above factors in the following way. 
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We use the population synthesis model GISSEL96, which is a revised version of Bruzual 
and Chariot (1993) (see also Leitherer et al. 1996), to calculate an SED from stellar components. 
The SED of the galaxy is then determined by its star-formation history. The SEDs of local 
galaxies are well reproduced by models whose star-formation rate declines exponentially (r 
model); i.e., SFR(t) oc exp(— t/r), where t is the age of the galaxy and r is the time scale of 
star formation. It is noted that different combinations between t and r generate a similar shape 
of SED in some cases. In this work, we therefore used r = 1 Gyr models with Salpeter's initial 
mass function (the power index of x = 1.35 and the stellar mass range of 0.01 < m/M < 125) 
to derive various SED types. We adopted the solar metallicity, Z = 0.02. For these models, we 
calculated SEDs with ages of t = 0.1, 0.5, 1, 2, 3, 4, and 8 Gyr, and called them SED1, SED2, 
SED3, SED4, SED5, SED6, and SED7, respectively; note that the SED templates derived by 
Coleman et al. (1980), elliptical galaxies (the bulges of M 31 and M 81), Sbc, Scd, and Irr 
correspond to those of SED7 (t = 8 Gyr), SED6 (t = 4 Gyr), SED5 (t = 3 Gyr), and SED4 (t = 2 
Gyr), respectively. We did not use an SED whose t was larger than the age of the universe 
at z. For our adopted model of the universe, the age of the universe exceeds 1, 2, 3, 4, and 8 
Gyr at z = 5.7, 3.2, 2.2, 1.6, and 0.6, respectively. For simplicity, we assumed that the dust 
extinction of the lensing galaxy to be Ay = 0. The effect of Ay on our discussion is mentioned 
below. As for the absorption by intergalactic neutral hydrogen, we used the average optical 
depth derived by Madau et al. (1996). In this procedure, we adopted an allowed redshift of 
between z = and those of quasars with a redshift bin of Az = 0.1. Figure 2 shows the 5-band 
absolute magnitude of galaxies corresponding to the observed limiting i'-band magnitude as a 
function of z and SED. 

3.2. Stellar Velocity Dispersion of the Lens Galaxy 

To estimate of the magnification factor by gravitational lensing, it is necessary to know 
the stellar velocity dispersion (a v ) of the lens galaxy. Given the redshift of a lens galaxy, we 
could estimate its absolute magnitude and then estimate a probable value of the stellar velocity 
dispersion. However, since there is no information about z, we tried to estimate the stellar 
velocity dispersion as a function of both the redshift and SED assuming that the Tully-Fisher 
relation established for a sample of galaxies in the local universe is also valid at a high redshift. 

Here, we used the Tully-Fisher relation derived by Sakai et al. (2000), 

M B = -8.07(logH/ 2 o- 2.5) -19.88, (1) 

where Mb is the absolute B magnitude and W20 is the full width at 20% of the maximum ve- 
locity. Assuming that a v = V vot /\^2 = H^o/^V^ where V rot is the rotation velocity, we estimated 
o~ v as a function of z and SED. We evaluated Mb using the method given in subsection 3.1. 
The results are shown in figure 3. In the case of SDSSp J103027. 10+052455.0, a v < 300 km 
s" 1 at z < 5.5. In the case of SDSSp J130608.26+035626.3, a v < 200 km s^ 1 at z < 5.5. 
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3.3. Magnification Factor by the Gravitational Lensing 

We are now ready to estimate the magnification factor by gravitational lensing. We 
adopt the singular isothermal sphere (SIS) model for simplicity (e.g., Binney, Merrifield 1998). 
In this model, the magnification factor for the brighter source can be expressed as 



where 9 is the angle between the lens galaxy and the source (i.e., SDSSp J103027. 10+052455.0 
and SDSSp J130608. 26+035626. 3 in this case), and #e is the Einstein angle, defined as 



where -Dls is the angular diameter distance between the lens galaxy and the source and -Dos 
is that between the observer and the source. M + of the brighter image and M_ of the counter 
image are related by M + /M_ = M + /(M + — 2). In figure 4, we show the magnification factor of 
the brighter source as a function of z and the SED types. Here, we adopt 9 = 1". The figures 
show that the maximum of M + is only M + ~ 5 (for SDSSp J103027.10+052455.0) and M + ~ 2 
(for SDSSp J130608. 26+035626. 3). On the other hand, M + is constrained from the condition 
that there is no counter image and results in M + < 2.31 (SDSSp J103027. 10+052455.0) and 
M + < 2.11 (SDSSp J130608. 26+035626. 3). The constraint that no lensing galaxy can be seen 
gives a similar result to that from the absence of any counter image. 

In the case of Ay > 0, the maximum M + , which is derived from the condition that there 
is no intervening galaxy, becomes larger. We emphasize, however, that the constraint on M + 
from the condition that there is no counter image is not changed. We therefore conclude that 
the two quasars are not strongly lensed. The gravitational lensing effect, if any, should be very 
modest. 

4. Concluding Remarks 

WL02 suggested that the SDSS z ~ 6 quasars are systematically brighter than expected. 
Namely, they investigated the probability that a quasar with a rest-frame 1450 A magnitude, 
M 1450 (AB), and a redshift of z was selected into the survey using the survey selection function 
given in Fan et al. (2001). They found that the cumulative probability of Mi45o(AB) for the 
z ~ 6 quasars is different from the expected one at the 95% significance level. However, they 
also suggested that this inconsistence is removed if the flux of one third of z ~ 6 quasars is 
overestimated by a factor of 10. 

In this paper, we cannot find any strong evidence that supports this suggestion. We have 
found no candidate of a lensing galaxy brighter than %' = 25.4 around the z=Q.28 quasar SDSSp 
J103027.10+052455.0 and i' = 25.8 around the z=5.99 quasar SDSSp J130608.26+035626.3. We 
have also found no counter image of the two quasars. We evaluated the possible magnification 
factor at the angle 9 = 1" using the SIS model: 1) M+ < 5 for SDSSp J103027. 10+052455.0, 




(3) 
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and 2) M+ < 2 for SDSSp J130608. 26+035626.3 in the case of A v = 0, and to be large in the 
case of Ay > 0. However, the constraints on M + from the condition that there is no counter 
image shows that M + is not very large in any case of Ay. Therefore, we conclude that both 
quasars are not strongly gravitationally magnified. 

Here, we mention a limitation on our analysis. If the angle between a quasar and a lensing 
galaxy (and a counter image) is much smaller than the seeing size, we cannot find evidence 
of gravitational lensing from our observations. High-resolution deep imaging at wavelengths 
shorter than the Lyman break of the quasar may give further constraints to gravitational lensing 
of this type. 

It is known that the lensing optical depth increases with increasing the redshift, and thus 
objects at higher redshift are more affected by gravitational lensing from a statistical point of 
view (e.g., Turner 1991; Barkana, Loeb 2000; WL02). Fan et al. (2003) noted that even without 
a detectable second image under high resolution, the possibility that luminous z ~ 6 quasars are 
modestly magnified by foreground galaxies is still high. As proposed by WL02, it will be very 
important to carefully search for lensing galaxies towards all high-z quasars by deeper imaging 
with higher resolution to constrain the magnification factor. To search for lensing galaxies 
that are overlapped with a quasar, high-resolution spectroscopy of the candidate of magnified 
quasars will also be important. 

We would like to thank the Subaru Telescope staff for their invaluable help. We also 
thank T. Hayashino for his invaluable help. This work was financially supported in part by the 
Ministry of Education, Culture, Sports, Science and Technology (Nos. 10044052 and 10304013). 
TN and MA are JSPS fellows. 
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Fig. 2. Distribution of the absolute magnitude of galaxies, which corresponds to the limiting magnitude, 
as a function of the SED and the rcdshift. 




Fig. 3. Distribution of the maximum velocity dispersion of the lens galaxy with > mum as a function of 
the SED and the rcdshift. 
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Fig. 4. Distribution of the maximum magnification factor by the lens galaxy with > m\i m as a function 
of the SED and the rcdshift. 



Table 1. Journal of imaging observations. 



Quasar 


Obs. date (UT) 


Total exp. (s) 


mum(AB)* 


^■quasar (AB) 


fwhmL (") 


SDSSp J103027.10+052455.0 


2002 April 14 


900 


25.4 


23.23±0.07 


1.2 


SDSSp J130608. 26+035626.3 


2002 April 14 


1200 


25.8 


21.95±0.03 


1.8 



* The limiting magnitude (3a) with a 2" aperture. 

t The full width at half maximum of stellar objects in the final image. 
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